Orbital-selective Mott transitions in the degenerate Hubbard model 
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We investigate the Mott transitions in two-band Hubbard models with different band widths. 
Applying dynamical mean field theory, we discuss the stability of itinerant quasi-particle states in 
each band. We demonstrate that separate Mott transitions occur at different Coulomb interaction 
strengths in general, which merge to a single transition only under special conditions. This kind of 
behavior may be relevant for the physics of the single-layer ruthenates, Ca2-iSrj;Ru04. 



Strongly correlated multi-orbital electron systems are 
among the most active topics in condensed matter 
physics. In the Mott insulators the addition of orbital 
to localized spin degrees of freedom leads to complex 
ordered phase diagrams. In itinerant electron systems 
multiple Fermi surface sheets appear with very distinct 
properties. Subtleties occur when localized and itinerant 
electrons coexist, as is well known in the case of itinerant 
d- and localized /-electrons which give rise to the rich 
physics of heavy Fermion compounds. In view of their 
very different character this is not so surprising. We may 
ask, however, whether it is possible to find coexistence 
of itinerant and localized electrons for degenerate non- 
hybridizing orbitals with small bandwidth difference. 

A nearly degenerate c?-electron system where multi- 
orbital properties obviously play an important role is 
the single-layer isovalent ruthenate alloy Ca2-a;Srj;Ru04. 
The end-member Sr2Ru04 is a well-known unconven- 
tional superconductor P, Q , while Ca2Ru04 is a Mott- 
insulating S = \ antiferromagnet [3|, The relevant 
4d-orbitals belong to the i2g-subshell. The planar struc- 
ture leads to very weak hybridization between orbital 
which have even {dxy) and odd parity {dyz,dzx) under 
the reflection z ^ —z. The complex evolution between 
these different end-members has led to various theoretical 
investigations 0, , and among others to the proposal 
that some of the d-orbitals display localized spin and 
orbital degrees of freedom, and others provide itinerant 
electrons. This orbital-selective Mott transition (OSMT) 
could explain the experimental observation of a localized 
spin S ~ 1/2 in the metallic system at x ~ 0.5 which 
is difficult to obtain from the entirely itinerant ab initio 
description 

The concept of the OSMT was recently challenged, in 
particular, by Liebsch whose dynamical mean field theory 
(DMFT) calculations suggested that two bands of differ- 
ent width coupled by electron-electron interactions would 
always undergo a common Mott transition 0- The aim 
of this paper is to revisit this problem and to analyze 
the Mott transition in the degenerate two-orbital Hub- 
bard model with different bandwidths. Using another 
form of DMFT, we find a different result and an OSMT. 
In addition, we show also that correlations can stabilize 
a commensurate filling of one band even when the total 



electron count is fractional. 

We consider the following Hubbard Hamiltonian with 
two orbitals. 
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where c[^„{ciaa) creates (annihilates) an electron with 
spin a{=] , I) and orbital index a{= 1, 2) at the ith site. 
We restrict ourselves to the case of non-hybridizing or- 
bitals, relevant to the ruthenates, and t-"^ denotes the 
hopping integral for orbital a, n the chemical potential, 
U [U') the intraband (interband) Coulomb interaction 
and J the Hund coupling. In the following, we restrict 
our discussions to the condition U = U' + 2 J, obtained 
by symmetry arguments for degenerate orbitals. 

We examine the stability of the metallic ground state 
of this model by means of DMFT which maps the lat- 
tice model to the problem of a single impurity con- 
nected dynamically to a "heat bath" The electron 
Green's function is obtained via the self-consistent so- 
lution of this impurity problem. We represent the two 
electron bands by semi-circular density of states (DOS), 
Pa{x) = 2/7rDa\/l — {x/DaY where 2D a is the band- 
width. For the case of identical hopping integrals for the 
two bands, t^^^ = tij (Di ~ D2), the role of orbital fluc- 
tuations has been discussed by means of DMFT [1E3- 
There are various methods to solve the effective impurity 
problem. We use here the exact diagonalization method 
proposed by Caffarel and Krauth lll| , since the quantum 
Monte Carlo (QMC) simulations used by Liebsch ^7] suf- 
fer from sign problems at low temperatures, in particular, 
if the Hund coupling is included. Additionally, we apply 
the hnearized version of DMFT (two-site DMFT) ^.ij, 
which allows us to discuss electronic properties well even 
in the vicinity of the critical point. We restrict our dis- 
cussions to the paramagnetic case to clarify the nature 
of the Mott transition. 

We first consider the case /xi = ^2 ^ U /2 + U' — J /2, 
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i.e. both bands are half filled. The quasiparticle weight 
Za, defined by = 1 — c?Re[I]ct(ti;)]/dw in terms of the 
self energy of each band, will be used to charac- 

terize the stability of the metallic state of the two bands. 
The results obtained with fixed ratios U' /U and J/U 
are shown in Fig. for half-filled bands. We first fo- 
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FIG. 1: The quasiparticle weights Zi and Z2 at half filling 
as a function of the Coulomb interaction U : (a) U' /U = 1.0 
(J = 0) and (b) U'/U = 0.5 {J/U = 0.25). The bandwidth 
is set as Di — 1.0 and D2 — 2.0. Open (closed) circles rep- 
resent the results for orbital a — 1(2) obtained by solving 
the DMFT impurity problem by means of the exact diagonal- 
ization of a small cluster (A'^ = 6). Solid triangles represent 
the Mott-transition points obtained by the two-site DMFT 
method, which produce the values quite consistent with those 
of the numerical diagonalization. Insets show the same plot 
for bandwidths Di = 1.0 and D2 — 5.0. 

cus on the case of C/ = C/' and J ~ [shown in Fig. 
(a)] with bandwidths Di = 1.0 and D2 = 2.0. When 
the Coulomb interaction is turned on, the quasiparticle 
weights Zi and Z2 decrease from unity in slightly differ- 
ent ways reflecting the difference of the bandwidth. A 
strong reduction of the quasiparticle weight appears ini- 
tially in the narrower band. However, when the system 
approaches the Mott transition, the quasiparticle weights 
merge again displaying a very similar dependence on U, 
and eventually reach zero at the same critical point. The 
inset shows the more extreme case of Di = 1.0 and 
D2 = 5.0 (very wide band). The common Mott tran- 
sition originates from the enlarged symmetry inherent 
in U — U' and J = 0, which will be discussed below. 
This result is in agreement with the conclusion by Lieb- 
sch 0. For small interaction strengths the decrease of 
the quasiparticle weight reflects the effective Coulomb 
interactions U /Da which are different for two bands of 



different width. Da- The dependence of the initial de- 
crease in Za on U/ Da reflects the general difference in 
Za- On the other hand, in the vicinity of the quantum 
critical point at the Mott transition the effect of the bare 
bandwidth is diminished due to the strong renormaliza- 
tion of the effective quasiparticle bandwidth allowing Zi 
and Z2 to vanish together [isl . 

The introduction of a finite Hund coupling J makes 
U ^ U' and leads to a qualitatively different behavior, as 
seen in Fig. (b). With increasing U keeping the ratio 
U'/U = 0.5 fixed, the quasi-particle weights decrease dif- 
ferently and vanish at different critical points: Ud ~ 2.6 
for Zi and Uc2 ~ 3.5 for Z2. Therefore, we observe an in- 
termediate phase with one orbital localized and the other 
itinerant, though strongly renormalized {Z2 ^ 1). The 
analogous behavior is observed for different choices of the 
bandwidths, if J takes a finite value [inset of Fig. ^(b)]. 
Although it is difficult to precisely determine the second 
critical point Uc2, this result certainly suggests the exis- 
tence of the OSMT with Uc2 > Ud- 

In Fig. 12 we show how the quasiparticle states evolve 
and then disappear inside of the Mott-Hubbard gap. The 
DOS is computed by the two-site DMFT scheme 0. 
In both cases the Mott-Hubbard gap develops as U in- 
creases and is accompanied by narrow quasiparticle mid- 
gap bands. For case (a) with J = 0, these quasiparti- 
cle bands disappear simultaneously, whereas for case (b) 
with finite J, they have different critical points, consis- 
tent with the results mentioned above. 

Repeating similar DMFT calculations for various 
choices of the parameters, we derive the ground-state 
phase diagram shown in Fig. 13 which displays some re- 
markable features. First, the metallic phase (I) remains 
stable up to surprisingly large Coulomb interaction U 
when U ^ U' (small J) . Here the Mott transitions merge 
to a single transition. This behavior originates from the 
high symmetry when U = U' ( J = 0) with six degener- 
ate two-electron onsite configurations: four spin config- 
urations with one electron in each orbital and two spin 
singlets with both electrons in one of the two orbitals. 
The additional symmetry in orbital/spin degrees of free- 
dom enlarges the phase space for charge fiuctuations and 
leads to a decrease of the Mott-Hubbard gap Eg at large 
U. A rough estimate of Eg can be obtained from the 
second moment of the hopping Hamiltonian for a state 
I a) with an extra electron (or hole). In (a|-ff^|a) = 
all possible configurations with the same onsite energies 
are considered as intermediate states 0. Because all 
charge excitations mix with each other, there is only one 
gap. Assuming a staggered spin or orbital configuration 
as the most dominant local correlation for neighboring 
two-electron sites we obtain for the effective hopping ma- 
trix element of an extra carrier T — y^f + tf . This al- 
lows us to estimate the Mott-Hubbard gap at large U: 

Eg = U - 2zT = U - 2z^tl + tl (2) 
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FIG. 2: The density of states pct{uj) at half fiUing: left (right) 
panel for a = 1 (2). (a) U' /U = 1.0 (J = 0) and (b) 
U'/U = 0.5 {J/U = 0.25). It is clearly seen that the Mott 
transitions occur simultaneously in (a), while the orbital se- 
lective transitions occur in (b). 
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FIG. 3: Phase diagram for two-orbital Hubbard model with 
Di — 1 and D2 ~ 2. In the phase (I) [phase (II)], both bands 
are in metallic (insulating) state. The phase (III) is induced 
by the orbital-selective Mott transition, where the metallic 
state coexists with the Mott insulating state. Since we are 
concerned with the ferromagnetic Hund coupling, J > 0, the 
relevant region in the diagram is U > U' . 



where z is the coordination number. We can get a simple 
estimate of Uc by setting Eg ^ leading to values Uc = 
2D1V2 (for Di = D2) and Uc = 2L>iV5 (for 2Di = D2 ). 
Both estimates are enhanced relative to the single-band 
case C7c = 2D [13. 

Away from the symmetric limit, i.e. U > U' (2J — 



U — U') orbital fluctuations are suppressed and the spin 
sector is reduced by the Hund coupling to three onsite 
spin triplet components as the lowest multiplet for two- 
electron sites. Applying the same scheme as above, we 
recognize that charge excitations with two electrons in 
one or the other of the orbitals do not mix, since all 
hopping processes included in preserve orbital con- 
figurations in the lowest multiplet sector. The effective 
hopping for each orbital is now Ta — ta/V^ assuming a 
staggered spin-1 state at half-filling. The reduction com- 
pared to the single band case occurs due to the locking of 
the spins into an onsite spin triplet. If we consider again 
the case 2Di — D2 we find two separate Mott transitions 
with critical values 



1 U' 

Uci = "^^1 + Y ^^'^ 



Uc2 = V2D1 + y . (3) 



In between the two transitions we find the metallic inter- 
mediate phase (III) with one band localized and Mott- 
insulating and one band itinerant. Within our DMFT 
scheme we have also confirmed that various choices of 
bandwidths give rise to the qualitatively same structure 
of the phase diagram as shown in Fig. |3| 

Our result for the Mott transitions is different from 
that of both Anisimov et al. and Liebsch 0, Q| . The for- 
mer group derived an OSMT for a special model which 
includes only the intraband Coulomb repulsion U within 
the DMFT approach, and drops effects due to coupling 
between the orbitals. This scheme was criticized by Lieb- 
sch, who took into account U, U' and J. He claimed 
based on a DMFT analysis that only a single Mott tran- 
sition occurs in the generic case. Liebsch's solution of 
the single-impurity problem within DMFT is based on 
QMC and iterative perturbation methods. The former 
suffers from sign problems which limit its validity at low 
temperatures, while the latter is an extrapolation from 
the small-C/ regime. Our DMFT analysis which uses the 
exact diagonalization of the impurity problem on a finite 
cluster is valid at zero-temperature and is not restricted 
to weak coupling. Our results show separate Mott tran- 
sitions occurring generically except for the special case 
with high symmetry U — U' (J = 0), for which the tran- 
sitions merge irrespective of the different bandwidths. 

We have so far treated the case of two individually 
half-filled bands. We now address the question what will 
happen when the electron count is non-stoichiometric. 
This problem may provide another key to understand 
the orbital selective Mott transitions in Ca2-xSrxRu04 
0, since each of the three original metallic bands pos- 
sesses fractional fillings. In order to study this kind of 
system, we introduce a finite hole doping, and observe 
how commensurability can emerge due to interactions. 

In Fig. 21 (a) we show the DOS which is computed by 
using the two-site DMFT. With increasing interactions 
quasiparticle states with large DOS appear around the 
Fermi energy in both bands. Enhancing the interactions 
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further we drive the first band insulating. This is in con- 
trast to the single band system, where finite hole doping 
obscures the Mott transition and always gives metallic 
behavior. In the two-band system, however, commen- 
surability in one of the bands gradually emerges, as is 
clearly seen in Fig. 0] (b). The electron number for 
the first band ni is plotted here. When U — 0, ni and 
n2(= I — S — Til) are smaller than 0.5 because of finite 
hole-doping (S = O.I). Coulomb interaction causes elec- 
tron transfer from one orbital to the other, giving rise 
to one half-filled band at a certain interaction strength, 
thereby causing an OSMT. 

(a) 1 
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FIG. 4; (a) The density of states pi{u>) and p2{ijj) for finite 
hole doping, 5 = 0.1. The Coulomb interaction is chosen as 
U = Q, 2.0, 3.0 and 4.0 (17' = 0.5f/ and J = 0.2517) from the 
top to the bottom, (b) The number of electrons in the orbital 
(ct = I) as a function of U when U' — 0.5U and J = 0.25?7. 



In conclusion we have discussed the Mott transitions 
in the degenerate Hubbard model with non-hybridizing 
orbitals which have different bandwidths by means of 
DMFT. Our analysis has shown that a single Mott tran- 
sition occurs when the Hund coupling is absent ([/ — U'), 
rendering the different bandwidth essentially irrelevant at 
the transition point, as discussed by Liebsch. In the more 
generic situation with finite Hund coupling, however, we 
find the OSMT. This is also true for non-stoichiometric 



systems. We believe that our study resolves the appar- 
ent contradictions on this issue. Moreover it sheds light 
on the nature of such Mott transitions, which may be 
relevant for the physics of the ruthenates. 
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